function c = vec_cross(a, b)
    c = [a(2)*b(3) - a(3)*b(2) ;
         -a(1)*b(3) + a(3)*b(1) ;
         a(1)*b(2) - a(2)*b(1)] ;
end